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Quasi-stationary criticality of the Order-Parameter of the d = 3 Random-Field Ising 
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The critical exponent f} — 0.17 ±0.01 for the three-dimensional random-field Ising model (RFIM) 
order parameter upon zero-field cooling (ZFC) has been determined using extinction-free magnetic 
x-ray scattering techniques for Feo.85Zno.15F2. This result is consistent with other exponents deter- 
mined for the RFIM in that Rushbrooke scaling is satisfied. Nevertheless, there is poor agreement 
with equilibrium computer simulations, and the ZFC results do not agree with field-cooling (FC) 
results. We present details of hysteresis in Bragg scattering amplitudes and line shapes that help 
elucidate the effects of thermal cycling in the RFIM, as realized in dilute antiferromagnets in an 
applied field. We show that the ZFC critical-like behavior is consistent with a second-order phase 
transitions, albeit quasi-stationary rather than truly equilibrium in nature, as evident from the large 
thermal hysteresis observed near the transition. 

PACS numbers: 61.10.Nz, 75.40.Cx, 75.50.Ee, 75.50.Lk 



INTRODUCTION 



One of the classic models of a phase transition in an 
intrinsically disordered system is the three-dimensional 
(d = 3) random-field Ising model (RFIM). After nearly 
three decades of progress in experiments, simulation and 
theory, the phase transition in this system is not yet 
well understood and remains a challenge to our funda- 
mental understanding of the statistical physics of dis- 
ordered systems. Extensive Monte Carlo simulations 
and exact ground state calculations have provided ev- 
idence for equilibrium critical exponents, including a 
value for the order-parameter exponent close to zero, 
but with a specific heat exponent that is not yet well 
determined 0, S S S ■ The specific heat exponent 
value, however, has been determined in experiments in 
the dilute antiferromagnet FexZni_xF2 at various mag- 
netic concentrations a; |^ 0| to be very close to zero 
upon warming through the phase transition. This is indi- 
cated by the apparent symmetric, logarithmic divergence 
as well as the field scaling of the optical birefringence, 
Faraday rotation, and specific heat amplitudes Q. The 
experimental determination of the order parameter has 
proven more elusive because the large number of vacan- 
cies at low magnetic concentrations allows the system 
to easily form domains that obscure the order-parameter 
critical behavior. 

An avenue for determining the order-parameter ex- 
ponent opened once it was understood that at higher 
magnetic concentrations these domains do not form. 
In the best studied system, the dilute antiferromagnet 
FexZni_xF2 in a field applied along the easy axis, do- 
mains do not form for magnetic concentrations x > 
0.75 0, 0, 0. As we will discuss below, the exper- 



imental results nevertheless do not agree well with the 
equilibrium Monte Carlo and exact ground state calcu- 
lations. Indeed, the experimentally-observed transition 
is clearly not in equilibrium, although critical behavior 
can be determined under conditions of monotonically in- 
creasing temperature. In this study, we describe the ap- 
parent critical behavior and characterize the hysteresis 
observed upon crossing the phase boundary. We also dis- 
cuss temperature reversals just below the phase bound- 
ary. In this way, we describe the unusual circumstance of 
self-consistent critical-like behavior in a clearly nonequi- 
librium system. The discrepancy with simulations is a 
consequence of the simulations being done under equi- 
librium conditions, something apparently not realized in 
the macroscopic experimental system. 

The phase transition in the d = 3 RFIM system 
Feo.85Zno.15F2 has been characterized in great detail for 
the zero-field-cooling (ZFC) procedure in which the sam- 
ple is cooled in zero field, the field is raised, and the sam- 
ple is warmed in constant field across the phase bound- 
ary. The transition appears to be second-order under 
ZFC; all of the critical exponents, measured at very small 
reduced temperatures, are self-consistent in that they ap- 
pear to satisfy the Rushbrooke equation, as described be- 
low. Nevertheless, it is well known that different behavior 
is observed upon field-cooling (FC), in which the sample 
is cooled across the phase boundary in the field. This 
has been observed for this sample in x-ray scattering |l2l | , 
neutron scattering optical birefringence 01 and op- 
tical Faraday rotation experiments. Hence, the ZFC 
phase transition is extraordinary in that it appears to 
be second-order if the temperature reversals are avoided 
but, in light of the hysteresis, it does not take place un- 
der equilibrium conditions as would be the usual case. 
To better characterize the apparent RFIM transition, the 
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FC critical behavior, and the hysteresis in general, needs 
to be explored more fully. We have done this for the 
order parameter in Feo.85Zno.15F2 using magnetic x-ray 
scattering. 

Characterization of the order parameter, the staggered 
magnetization (Ms), is achieved by determining the tem- 
perature dependence of the antiferromagnetic Bragg scat- 
tering intensity, which is proportional to AIs, versus the 
temperature T. The order parameter is expected to be- 
have as 



Ms = Mot 



(1) 



for i << 1, where t = {Tc{H) - T)/Tc{H) is the reduced 
temperature and Tc{H) is the transition temperature. As 
discussed in a previous Letter , which focused on the 
ZFC measurements, neutron scattering techniques can 
not be used to characterize the order parameter criti- 
cal behavior in high-crystalline-quality bulk crystals be- 
cause of severe extinction effects, which tend to modify 
the temperature dependence of the Bragg intensity as 
the temperature is lowered. The effect arises when the 
scattering sample region is so thick and the crystal is so 
perfect that the beam is selectively depleted of neutrons 
that satisfy the Bragg condition for scattering. This pre- 
vents an accurate determination of the exponent (3. The 
x-ray technique, on the other hand, is essentially free of 
extinction effects and the order-parameter criticality can 
be measured accurately near Tc{H) 0|. This is possi- 
ble since the thickness of the scattering region is limited 
by the strong temperature-independent charge scatter- 
ing and, within this region, the beam is never depleted 
of x-rays meeting the Bragg condition since the magnetic 
scattering is weak. 

A magnetic concentration near the one in this study, 
X = 0.85, is crucial to this order parameter characteriza- 
tion. Many prior attempts 0, ^3 to determine the crit- 
ical behavior of the order parameter proved unsuccessful 
because x < Xv, where x^ = 0.754 is the magnetic va- 
cancy percolation |inl |. Domain formation obscures the 
RFIM critical behavior below the transition at Tc{H) in 
FexZni_xF2 and its less anisotropic (i.e., smaller Ising 
anisotropy) isomorph MnxZn2_xF2 for x < x^. The con- 
centration X — 0.85, while being greater than x^, is low 
enough to generate significant random-field effects at the 
maximum field of our experiment, H = 11 T. It is, at 
the same time, high enough to avoid the complication of 
contributions to the scattering intensity from fractal-like 
percolating vacancy structures that appear close to 



EXPERIMENTAL DETAILS 

We are interested in the magnetic scattering intensity 
near the antiferromagnetic zone center (10 0). Hence, we 
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FIG. 1: The energy of the x-ray was carefuUy fine-tuned to 
ensure that the effect of multiple scattering is minimized near 
the (10 0) antiferromagnetic Bragg point, which occurs near 
9 — 5.1 degrees. The temperature is several degrees above 
the transition temperature, so no magnetic peak is present. 
The multiple scattering peaks appear symmetrically on either 
side of the Bragg point. These data demonstrate that the 
shoulder peaks move apart and become weaker as the energy 
is tuned close to 13605 eV. The curves are Lorentzian fits. 
The widths increase as the peaks move apart and decrease in 
amplitude, allowing the Bragg peak to be clearly discernible 
at low temperatures, as shown in Fig. 2. 



need not worry about strong charge scattering contribu- 
tions and can easily discern the relatively weak magnetic 
response. The measurements were made at the high- 
field magnet facility on beam line 7-2 of the Stanford 
Synchrotron Radiation Laboratory. A monochromatic 
x-ray beam was obtained using a Si(lll) double-crystal 
monochromator from a spectrum produced by a wiggler 
insertion device. The x-ray energy was carefully tuned 
to a value between 13.5 and 14 keV to minimize the ef- 
fect of energy-sensitive multiple scattering peaks around 
the magnetic Bragg point 'XM- T^^^ incident x-ray en- 
ergy was well-defined, within about 10 eV, whereas the 
detector had a half-width energy resolution of 300 eV. 
As shown in Fig. ^ at temperatures above the transition 
temperature, where no magnetic scattering occurs, the 
x-ray energy was adjusted until the multiple scattering 
peaks move apart as far as possible ^1^- From Fig.[21it is 
clear that the multiple scattering peaks do not affect the 
analysis of the magnetic peak below the transition for an 
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FIG. 2: Typical scattering at low temperature for H — 0. The 
energy is adjusted to minimize the multiple scattering, which 
is nevertheless still present, but well separated from the Bragg 
peak. The magnetic Bragg peak is about half as wide as the 
multiscattering peaks, as determined from the Lorentzian fits 
shown in the figure. 



x-ray energy of 13595 eV. 

The sample has a finely poUshed face, approximately 
13 mm^ in area, and is 0.9 mm thick, with the a-axis 
perpendicular to the poHshed face. It was mounted such 
that the c-axis was along the vertical field. The well 
polished surface prevented spurious effects such as those 
observed in previous experiments 0| , where a transition- 
like behavior was observed that disappeared when that 
the faces were polished. Presumably, this is a result of 
strong pinning at the locations of the scratches which 
prevented RFIM correlations to develop. The chemical 
homogeneity of the crystal was determined by a room 
temperature optical linear birefringence technique to 
be 0.45%/cm. The rounding of the transition is approx- 
imately 50 mK, corresponding to a tiny reduced tem- 
perature of 0.0004 on either side of the transition. The 
crystal was mounted on a thin silicon bar of dimensions 
0.8x 1.5x 15mm'^ and placed in an atmosphere of approxi- 
mately 10-20 mbar of Helium gas to achieve stabilization 
of the sample temperature to within approximately 10 
mK. The zero-field transition temperature was measured 
to be Tn = 66.7 K, consistent with birefringence mea- 
surements on the same sample 7J smd with a magnetic 
concentration x = 0.85 17]. 

The lattice parameters of the sample were determined 



to be approximately a — 4.68 A and c = 3.27 A 
near the transition temperature. The fuU-widths-at-half- 
maximum measured at the (1 0) magnetic Bragg point 
were 4 x 10~^, 4 x 10"'^, and 4 x 10"'^ reciprocal lattice 
units (r.l.u.) for the transverse, longitudinal and verti- 
cal directions, respectively. Three conventional thermal- 
cycling procedures were employed: ZFC; FC; and field- 
heating (FH). As described previously, in ZFC the sam- 
ple is cooled across the phase boundary with H = and 
data are taken while warming with H > 0, whereas in 
FC the data are taken while cooling across the phase 
boundary with H > 0. In the FH procedure, the sam- 
ple is first cooled through the transition in a field, and 
then the data are taken at that field value at successively 
higher temperatures. During each procedure, when the 
sample temperature was not continuously changed, the 
sample was held at each temperature for at least 20 min 
before taking data to ensure that the temperature and 
system stabilized. The temperature was stabilized be- 
fore q scans were obtained. However, when monitoring 
the peak intensity, it was more convenient to let the tem- 
perature change continuously. The data taken at H = 
and iJ = 11 T in this study, except where noted, were 
obtained in transverse (I q 0) scans typically consisting 
of 41 points, about 15 of which covered the Bragg peak. 
At each point, the intensity was counted for 30 to 45 sec- 
onds, depending on the temperature of the scan. At other 
fields, Bragg intensities were obtained at g = only. 

The antiferromagnetic transition for x = 0.85 has 
been shown to be stable in applied fields as high as 
H = 18 T To make sure that the ordered sys- 

tem is well-behaved at low T as we raised the field, we 
first cooled the sample in zero field to 20 K. The mag- 
netic field was then slowly raised to 11 T at a rate of 0.4 
T/min for < 9 T and 0.1 T/min for 9 < i? < 11 T. 
The peak intensity of the order parameter was monitored 
as we raised the field. Another set of measurements was 
similarly taken at T = 45 K. The results at both temper- 
atures are shown in Fig. O The intensity at T = 20 K 
remains essentially unchanged as the field is increased 
since the data are all taken deep within the ordered re- 
gion. The result for 45 K shows a slight field dependence 
since, for this case, the phase boundary (at 63.7 K for 
= 11 T) is approached more closely, although it is 
still quite far away. Those data demonstrate the stabil- 
ity of the order, i.e., no apparent phase boundaries were 
crossed. No change in the line width of the Bragg scat- 
tering was observed upon application of the field, further 
attesting to the stability of the order at low temperatures. 

The x-ray scattering technique has the advantage of 
very high momentum resolution, allowing a detailed 
study of the Bragg peak line shapes. It is important to 
use fine collimating slit widths to observe details of the 
scattering line shapes. For this purpose, the horizontal 
slits were configured to be 0.5 mm, approximately 30 cm 
in front of the sample and 1.1 mm the same distance after 
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FIG. 3: Field dependence of the Bragg peak intensity mea- 
sured at T = 20 K and 45 K after cooling in ii" = 0. The 
long-range order remains stable in fields up to = 11 T, i.e., 
no apparent phase boundaries are crossed. 

the sample. At approximately 70 cm behind the sample 
is a second slit which was set at 1.35 mm. Care was taken 
to ensure that the sample was well aligned for all mea- 
surements while using this narrow slit configuration. For 
measurements of the peak intensity versus temperature, 
the slits after the sample were opened slightly to reduce 
the sensitivity of the intensity to the precise alignment 
of the sample and beam so that wide ranges of T and H 
could be accessed without continual adjustments. 

Figure 0] shows typical transverse scans at different 
temperatures in a field H = 11 T after ZFC. The 
peak line shapes appear to be consistent with resolution- 
limited Gaussians for all temperatures below the tran- 
sition temperature for H = 1\ — 63.7 K. As the 
temperature approaches Tc, the intensity gradually di- 
minishes. The relatively very weak critical scattering is 
not apparent in this figure. 

QUASI-STATIONARY CRITICAL BEHAVIOR 

Figure |31 shows the (1 0) Bragg intensity at = 
and, for ZFC, dX H — \1T versus temperature, with the 
momentum and temperature independent background 
subtracted. The background depends on the precise ex- 
perimental configuration, but not on the thermal cycling 
procedure used in collecting data. The background is 



FIG. 4: Representative transverse scans for the different tem- 
peratures below Tc{H) — 63.7 K taken with H = 11 T after 
cooling in H = 0. Each scan is displaced vertically by 0.1 
units from the scan below it for clarity. The solid curves are 
results of least-squares fits to a Gaussian line shapes with a 
half-width-at-half-maximum equal to 2.1 x 10^* r.l.u. 

mostly from sources other than the crystal itself. For 
comparison of the background to the Bragg intensity sig- 
nal, typical background counts for the _ff = 11 T scans 
were eight counts per second whereas the q = inten- 
sity was 300 counts per second at T = 37 K. Above the 
transition, the background-subtracted intensity at small 
q results only from the critical scattering and goes to zero 
well above (H) , indicating that there are no discernible 
contributions from multiple scattering to the measured 
Bragg intensities. To determine the critical scattering 
for the H = and 11 T scans, neutron scattering line 
shapes, obtained with a sample of nearly the same mag- 
netic concentration |l3l | using a procedure described pre- 
viously [l^ Hij . were folded with the x-ray resolution, 
and the overall q = amplitude was adjusted to fit the 
H ^ 11 T data above Tc{H). Insets a and b in Fig. [S] 
show the critical scattering contributions for H ^ and 
11 T, respectively. As a result of the high momentum 
resolution of the x-ray technique, the critical scattering 
contributions at H = 0, which are nearly Lorentzian for 
H — iS, i3j almost negligible (inset a). For H > 
(inset b), however, the critical scattering has been shown 
to be non-Lorentzian and to have a much larger inten- 
sity at small q 19]. Recently, it has been shown that 
this scattering is consistent with fractal spanning clus- 
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FIG. 5: The Bragg ZFC intensity, Ib, versus T for a) H = 
and b) = 11 T with the momentum and temperature in- 
dependent background subtracted. The results in the insets 
show the critical scattering contributions to the x-ray inten- 
sity for H = and _ff = 11 T, respectively, determined from 
neutron scattering measurements, as described in the text. 



FIG. 6: The same ZFC data as in Fig. corrected for the 
critical scattering contribution, plotted as the logarithm of 
the intensity versus the logarithm of t. The solid line for 
H = 11 T indicates RFIM behavior with (3 = 0.17, while the 
solid line for H = reflects conventional random-exchange 
behavior with (3 — 0.35. 



ters, nucleated at random-field pinning sites, that form 
and grow as as the transition temperature is approached 
from above |l3| . Consequently, a small contribution to 
the (7 = scattering is more discernible for the 7? = 11 T 
data. After the critical scattering contribution is sub- 
tracted from the overall Bragg scattering intensity, the 
order parameter exponent can be determined using Eq. 
^ As shown previously |l2j |. the exponent determina- 
tion is not very sensitive to the details of the background 
subtraction. 

Although neutron scattering measurements using 
Feo.85Zno.15F2 ^ and Feo.93Zno.07F2 [H show no evi- 
dence for micro-domain structure formation in the crit- 
ical scattering, H > hysteresis in Ib is evident. FC 
intensities are larger than the ZFC ones, which is a re- 
sult of extinction. For ZFC samples, the ordering is so 
perfect that relatively few neutrons satisfy the Bragg con- 
dition. With FC, even with no large-scale domain struc- 
ture, there is apparently enough disorder to strain the 
crystal, through magnetostrictive effects, allowing more 
neutrons to scatter. The x-ray Bragg scattering also 
shows hysteresis, but in this extinction-free case the ZFC 
data are higher in intensity. FH data are intermediate 
between the ZFC and FC curves. We note that specific 



heat critical behavior measurements also show hysteresis 
very close to Tc (H) at this concentration P, Q . 

The Bragg intensity curves for > in Fig. O clearly 
approach Tc{H) vertically. This is characteristic of ex- 
periments 9, .13, .19] and simulations 10^ 20.1 fo r x > Xv 
and is in stark contrast with experiments |l4L lisl Elf 
and simulations for a; < Xv, where Ib approaches 
Tc{H) horizontally. The latter behavior is attributable 
to micro-domain structure formation, which is energeti- 
cally favorable when the vacancies percolate through the 
crystal, as shown in Monte Carlo simulations 10]. 

Figure IHl shows the logarithm of /b, with the con- 
stant background and critical scattering contributions 
subtracted, vs. the logarithm of t for H — and for 11 T 
under ZFC cycling. The values of Tc{H) were determined 
from fits to the data. For 0.0007 < t < 0.03 and H ^ 0, 
we find /3 = 0.35 ± 0.02 (lower solid line), which agrees 
well with several experimental and theoretical determi- 
nations for the random-exchange Ising model, as recently 
discussed ,21j,,2^,2^]. For 11 T, a crossover from random- 
exchange to RFIM critical behavior occurs near t = 0.03, 
consistent with birefringence measurements and the 
data can be fit to a single power law only in the range 
0.0001 < t < 0.03. The fit over this range yields the 
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FIG. 7: The logarithm of Bragg peak intensities Ib vs. tem- 
perature in fields of 7, 8.5, 10, and 11 T. All data have been 
corrected for the critical scattering contribution. Solid lines 
with a slope of 2/3 = 0.34 are added to the data set at each 
field. The amplitudes of the data sets are adjusted to agree 
at large t. 

exponent /3 = 0.17 ± 0.01 for H = 11 T, as indicated by 
the upper solid line in Fig. |S| 

In addition to the H = 11 T data described above, 
we also measured the order parameter at H ^ 7, 8.5, 
and 10 T to investigate crossover effects. All of the 
data were taken upon implementing the ZFC protocol. 
After we realigned the (1 0) Bragg peak position at 
iJ = 11 T using transverse scans, the sample was slowly 
warmed through the transition while the intensity, /b, 
was recorded. The rate of temperature change was con- 
trolled to be 0.2 K/min away from the transition, and was 
decreased to 0.02 K/min close to the transition. After the 
sample was warmed above the phase boundary and sub- 
sequently cooled, we re-examined the beam alignment to 
ensure that the peak intensity was measured. We found 
the peak position to change very little upon cycling. The 
same procedures were repeated for fields of 10, 8.5, and 
7 T. Figured summarizes the ZFC order parameter mea- 
surements at various fields. Three common features can 
be discerned. First, far below the transition tempera- 
ture, the Bragg peak intensity tends toward saturation 
since all magnetic moments in the beam-illuminated re- 
gion achieve long-range order at such fairly low temper- 
atures. Second, for reduced temperatures in the range 
< t < 10^°-^, the behavior of \ogiQ{I) is very sim- 



ilar to the REIM behavior, as shown in Fig. the dif- 
ference in the slope results from the variation of Tc with 
field, since the transition temperature is determined by 
the random-field behavior. Third, as the temperature 
approaches Tc{H), the logarithm of Ib in different fields 
exhibits the same exponent which is indicated by the par- 
allel straight lines through the experimental data. The 
fitted exponents of the order parameter at iJ = 7, 8.5, 
10 and 11 T are /3 = 0.18 ± 0.02, 0.18 ± 0.02, 0.17 ± 0.02 
and 0.17 ± 0.01, respectively, where the quoted errors 
are statistical. The higher fields provide larger ranges of 
reduced temperatures over which to fit the data and so 
yield more reliable values for p. From these results, we 
obtain a more precise order-parameter critical exponent, 
P = 0.17±0.01, than reported in our previous study [l^ . 

In Fig.|Slis plotted the logarithm of H versus the loga- 
rithm of t where REIM crosses over to RFIM critical be- 
havior. This value is defined, for the data shown in Fig. 
13 as the intersection of a straight line through asymp- 
totic RFIM behavior and one through the REIM behav- 
ior. From the scaling variable tH^^/'^ , we see that the 
slope in Fig. |S| yields the random- exchange to random- 
field crossover exponent (j). Indeed, the value obtained, 
(j) = 1.4±0.05, agrees well with the earlier experimentally 
determined and the most recent theoretically established 
values, both of which are 4) = 1.42 ± 0.02 for three di- 
mensions im m m mig . 

Through the Rushbrooke scaling relation 

2/3 7 a > 2 , (2) 

which is usually satisfied as an equality, [3 is related to the 
universal critical exponents a (for the specific heat) and 7 
(for the staggered susceptibility) of the d = 3 RFIM. The 
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FIG. 9: The temperature dependence of the Bragg peak in- 
tensities measured in the temperature range of 20 K and 60 K 
with H = 11 T. The dashed hne is a guide to the eye. 



experimentally determined specific heat peak is nearly 
logarithmic and very symmetric close to Tc{H), consis- 
tent with a = O.OOibO.Ol 0,0- Recent neutron scattering 
analyses 0, 0| yield a value of 7 = 1.68 ± 0.03, con- 
sistent with earlier results. Therefore, the experimental 
value (3 ~ 0.17 ±0.01 is consistent with Rushbrooke scal- 
ing 2f3 + J + a — 2.02 ± 0.06, even though the system is 
clearly not in equilibrium. 

None of the ZFC results are dependent on the rate 
at which the sample was warmed, typically between 
0.2 K/min and 0.02 K/min. However, the results are 
quite sensitive to temperature reversals of even a few mk, 
including overshoots of the set point when stabilizing the 
temperature. Such overshoots were meticulously avoided 
in the critical behavior measurements. 



THERMAL HYSTERESIS 

We next turn our attention to the effects of hysteresis 
and temperature reversals, including the difference be- 
tween the critical behavior observed upon ZFC and PC. 
The order parameter measurements display significant ir- 
reversibility in both T and H cycling procedures. We will 
address field hysteresis effects in the next section. 

A ZFC temperature cycle, followed by repeated cool- 
ing and heating at i? = 11 T, is shown in Fig. |51 for 
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FIG. 10: Behavior of the Bragg scattering upon tempera- 
ture reversals after ZFC for H = in (top panel), and after 
the sample was FC to 40 K (bottom panel). In the upper 
panel, ZFC Bragg peak intensities at = 11 T, as well as 
after several reversals in temperature below Tc = 63.7 K, are 
shown. The data were obtained in pairs of warming and cool- 
ing scans (w-c), except for the final set, which was obtained 
only upon warming (w). The curves are guides to the eye. In 
the bottom panel, intensities from the final warming through 
Tc shown (curve #6) in the previous figure are plotted along 
with FC intensities and reversal intensities a.t H = 11 T. The 
data shown were obtained in warming (w) or cooling (c) scans. 
The curves are guides to the eye. 



the temperature range between 20 K and 60 K. There is 
no observable hysteresis in this temperature range; the 
data taken upon cooling are essentially identical to those 
obtained upon heating. 

However, if we heat the sample to a temperature suffi- 
ciently close to the phase boundary, then, upon reversing 
the temperature, the intensity of Bragg scattering shows 
significant irreversibility, as shown in the upper panel of 
Fig. ^1 The data obtained upon cooling exhibit lower 
intensities compared to the data obtained upon the ini- 
tial ZFC warming. A significant discrepancy is observed 
after the sample is cooled to T = 55 K. The behavior 
of temperature reversal amplitudes of the order param- 
eter after cooling with H = and raising the field was 
further investigated by repeatedly cooling and warming 
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the sample, using rates described earlier. The top panel 
shows the measurements taken upon cycling after cooling 
with H = and raising the field without exceeding Tc 
at = f f T. All the reversals close to Tc take place at 
T = 63.40 K, which is 300 mK below T^. For reversal 
measurements without crossing the phase boundary, the 
amplitudes of the Bragg peaks gradually decrease after 
each cycle, but the warming curve follows the previous 
cooling curve quite well. In this way, the overall ampli- 
tudes slowly decrease upon cycling. After several cooling 
and warming loops, the amplitude is finally stabilized 
at the lowest cooling curve (scan ^^5) shown in the top 
panel. The last set of data (scan #6) was obtained by 
warming through the phase boundary. 

In the lower panel of Fig. ^1 the data obtained while 
warming across the phase boundary are repeated to serve 
as a reference. The sample was cooled using the FC 
procedure to 55.00 K. The same temperature loop mea- 
surements were then repeatedly taken between 55.00 K 
and 63.40 K. Contrary to the cycling shown in the up- 
per panel, the cooling and warming curves do not exhibit 
as clearly a slow approach to the stabilized equilibrium 
state from above, but instead seem to reach the stabi- 
lized behavior quickly from below upon warming. Fur- 
ther thermal cycling does not change the intensity. The 
initial FC curve is quite distinct from all other heating 
and cooling curves below Tc. The FC transition appears 
to be at the same temperature as the ZFC one, within 
experimental accuracy. 

In Figure ^2 the logarithm of the peak amplitude 
versus logarithm of t is plotted for various heating and 
cooling procedures. Only data of some typical scans are 
shown. Clearly, ZFC and FH curves for = 11 T have 
the same critical power- law behavior with P = 0.17 at 
small t. For temperatures outside the random-field criti- 
cal region, the data cross over from the RFIM asymptotic 
critical behavior to the random- exchange Ising behavior 
as t increases. The order parameter measured using the 
FC procedure shows strikingly different behavior from 
that of ZFC and FH ones; the overall intensities stay 
much lower and there is no observable crossover from 
REIM to the RFIM behavior. The exponent for the FC 
has a value of /5 « 0.35, but perhaps shows some round- 
ing at t smaller than 0.002. Over a large range of re- 
duced temperature, the results obtained upon FC seem 
very similar to the random-exchange Ising model. Two 
different cooling rates, 0.002 and 0.23 K per minute, were 
used for the FC protocol. Even though these rates differ 
by two orders of magnitude, we find that the measured 
intensity is essentially insensitive to the choice of cooling 
rate. In order to compare the H — and FC data, the 
zero-field peak intensity is multiplied by 0.70, which re- 
flects the difference of line shapes for the two cases, as 
we will discuss below. 

Thermal cycling loops at = 10 T after FC are shown 
in Fig. 1121 As the reversal temperature is lowered away 




logio(t) 



FIG. 11: A subset of the same data shown in Fig. 1101 plotted 
as the logarithm of the amplitudes versus the logarithm of t. 
Also included are data from a scan at _ff = as well as a 
fast FC scan. The fast scan was obtained by cooling through 
the transition at a rate of 0.23 K per minute. The data from 
scan 4^=1 are only those obtained upon warming. The top two 
curves represent ZFC and the next two are FC and heating 
after FC. They all indicate a critical exponent consistent with 
/3 = 0.17 ± 0.01. The solid line at the bottom has a slope 
indicating /? = 0.35. 



from Tc(-ff), the warming curves, above the flrst few 
points after the reversal, gradually display normal FH 
behavior. For data taken with the lower temperature re- 
versal point very close to Tc, there is little discernible 
difference between warming and cooling data. 

The irreversibilities of the random-field Ising model or- 
der parameter are also reflected in the peak intensity dif- 
ference between ZFC and reversal curves at low tempera- 
ture. The size of the intensity difference depends on how 
close one chooses to reverse the temperature near the 
phase boundary, as shown in Fig. 1131 Four points close 
to Tc = 62.85 K at 77 = 10 T were chosen to reverse the 
temperature: 62.10 K; 62.35 K; 62.55 K; and 62.62 K. 
Each reversal measurement was taken after ZFC prepa- 
ration. The intensities for ZFC curves are normalized at 
the lowest temperature so that they all have the same 
amplitude for comparison. For reversals close to Tc{H), 
the shapes of the reversal curves in the critical region are 
quite different from the ZFC ones. The shapes for the 
cooling curves suggest a larger /3 since they appear not 
to be as steep as the ZFC ones. 
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FIG. 12: Cycling measurements at = 10 T after FC 
with low-temperature reversals at four different temperatures, 
62.35, 62.20, 61.70 and 60.30 K. The warming curves gradu- 
ally exhibit normal FH features as the reversal temperature 
moves away from Tc, but the cooling curve is always lower. 
The triangles represent data under FC conditions and the 
crosses represent data under warming conditions. 

The temperature dependences of the x-ray scatter- 
ing line shapes are shown in Fig. 1141 under zero field 
cooled, field cooled and field heated conditions for H = 
10 T. Under ZFC, the data have typical Gaussian shapes 
with a half- width- at-half-maximum (HWHM) of = 2.1 x 
10~^ r.l.u. For FC and FH conditions, the line shapes and 
intensities are distinctly different from ZFC. The tails in 
the former cases are much larger and the central inten- 
sities are smaller. The temperature dependence of the 
widths for ZFC, FC and FH is shown in Fig. UHl All 
data were fit using Gaussian line shapes. The HWHM 
from Gaussian fits for the FC and FH line shapes is much 
larger than that for ZFC. Gaussian fits describe the data 
well and indicate 3.5 x lO"** r.l.u. for the widths under 
FC and FH. Although it is not known if Gaussian fits 
are the correct ones to use for FC and FH, they do work 
well and afford direct comparisons of widths for the three 
cases. 

In Fig. the peak intensities are plotted versus re- 
duced temperature. For a more meaningful comparison, 
we should compare the integrated intensity versus t. Fig- 
ure El shows the shape of transverse scans for ZFC and 
FH. In both cases, the line shapes appear to be Gaus- 
sian, but with different peak widths. If we take that peak 



FIG. 13: A series of reversal curves near the phase boundary 
at i/ = 10 T. The reversal points are close to, but below the 
phase boundary. 

broadening into account, the integrated intensity should 
be the product of peak width and peak intensity. Figure 
1 161 shows the logarithm of integrated intensity versus the 
logarithm of t. Indeed, within the experimental accuracy, 
the FH data collapse onto the ZFC data curves. 

Figure El shows the difference between FH and FC 
data sets increasing with the strength of the applied field. 
The hysteresis is difficult to discern in these measure- 
ments ioT H < 8 T, whereas for H = 10 T the different 
shapes for FH and FC are quite evident. 



FIELD HYSTERESIS 

Finally, field-cycle measurements were carried out to 
study the history dependence near T = 63.10 K, the 
transition temperature for H = 9 T. The sample was 
first warmed into the paramagnetic state at 70 K, then 
cooled in zero field to 63.10 K. The intensity was recorded 
as the field was slowly raised to 11 T, where Tc = 63.7 K, 
followed by a series of field cycling procedures with the 
temperature held at 63.10 K. The result is presented in 
Fig. El The general features of the experimental data 
are similar to those of the thermal cycle measurements 
in Fig. El The initial "field raising" data, prepared by 
cooling in zero field and represented by open triangles, 
stay the highest and the intensity goes to zero rapidly as 
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FIG. 14: Representative transverse x-ray scattering scans for 
the FC and FH protocols. Data were taken at // = 10 T. The 
solid curves are the results of least squares fits to a Gaussian 
shape. A scan taken under ZFC conditions is also plotted in 
the upper panel for comparison. 

Tc{H) is approached. As the field is lowered, significantly 
less of the peak Bragg intensity is recovered, as shown for 
data represented by the solid symbols. When the field is 
raised after first lowering it, the Bragg intensity initially 
has an intensity similar to that at low field, but eventually 
crosses over to the initial "field raising" curve, in the 
vicinity of He- The curvature is more like the initial "field 
raising" procedure near the transition. Similar behavior 
is seen upon raising the field beginning at H = 7.6 T. 

CONCLUSION 

In summary, we used extinction- free x-ray scattering to 
study the hysteresis of the order parameter of random- 
field Ising model. The critical exponent /3 = 0.17 ± 0.01 
is obtained for ZFC, whereas a different critical- like be- 
havior, very similar to the random-exchange behavior 
seen for H = 0, takes place for FC. The crossover 
from random-field Ising model to random-exchange Ising 
model was observed and the value <j> = 1-40 ± 0.05 ob- 
tained for the crossover exponent is consistent with the 
earlier studies and theory. We examined the history- 
dependent critical behavior in detail for various thermal 
and field cycling. It is clear from the different cycling 
experiments that the system under FC is not in equilib- 



FIG. 15: The HWHM of x-ray scattering line shapes for ZFC, 
FG and FH. ZFG line shapes are much narrower than those 
for FG and FH. In the graph, representative error bars are 
shown. 

rium. However, the critical behavior upon ZFC is rate 
independent and is consistent with a second-order tran- 
sition if the temperature is never reversed. The ZFC 
order-parameter critical behavior, i.e. the power law be- 
havior with j3 = 0.17 ± 0.01 and the crossover to REIM 
behavior, must be associated with rather stable, quasi- 
stationary states of the magnetic long-range component 
of the order. FH data exhibit line shape widths larger 
than ZFC and similar to FC, suggesting that some disor- 
der introduced upon FC remains. Aside from the slightly 
wider line shapes, the critical behavior of the FH data 
closely resembles the ZFC behavior. This suggests that 
it is the heating itself that is important in the manifes- 
tation of the new RFIM critical behavior and not the 
history of the temperature-field cycling. 

Although these findings are not in good agreement 
with equilibrium simulations and ground state calcula- 
tions, as mentioned earlier, they are in good agreement 
with non-equilibrium Monte Carlo studies \2^. Further- 
more, they are consistent with studies on uniaxial re- 
laxor ferroelectrics, a rather different experimental real- 
ization of the random-field Ising model |30l |31| in which 
the order-parameter critical exponent is observed to in- 
crease from (3 — 0.13 to [3 = 0.30, a value close to the 
random-exchange value, when the initial polarization is 
varied from 100% to 0.8%. It is argued for the ferroelec- 
tric system that a lower initial polarization corresponds 
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FIG. 16: The logarithm of the integrated scattering intensity 
of transverse scans vs. the logarithm of t for ZFC and FH 
protocols. The two sets of ZFC data represent experiments 
that were separated by several months in time. 



FIG. 17: Field dependence of the Bragg scattering for < 
H < lOT. The difference between FC and FH increases as the 
field increases. The FC intensity at low temperature becomes 
lower than the FH intensity. 



to a compensation of the random- field by domain walls. 
In the present case of the dilute antiferromagnet in a 
uniform field, it is clear that the large suppression of 
the transition is still present upon FC and that random 
fields are therefore not suppressed. However, it is also 
clear that the new critical behavior that is observed upon 
ZFC does not occur when the sample is FC. Mean-field 
treatments [s^, H^l , exact ground state calculations [s^l 
and Monte Carlo simulations indicate a very com- 
plex energy landscape and unusual characteristics near 
the RFIM phase transition. Perhaps in FC the fractal 
spanning cluster structures formed in our sample as it 
is cooled towards Tc{H) influence the character of the 
ordering process upon FC 0, ■ The spanning clus- 
ters form as the transition is approached from above and 
represent an ordering process quite different from pure 
systems. If the system cannot readily evolve from that 
configuration just above Tc{H) to long-range order be- 
low, this might indeed result in the severe hysteresis we 
have observed. In the ZFC process, this high magnetic 
concentration sample starts from a fully ordered lattice, 
retaining a single domain structure, and is perhaps less 
influenced by spanning cluster structures that form above 
Tc{H). Below Tc{H), metastability upon cooling has 
been described recently in terms of instantons which are 
a result of the complicated energy landscape due to the 



random fields [36j- Reversals of the temperature below 
the transition surely represent states in between those 
formed under ZFC and FC. A theoretical understand- 
ing of these differences and why they occur is lacking at 
this time. However, it is clear that this transition, with 
its coexistence of second-order-like critical behavior mea- 
sured to very small reduced temperatures as well as se- 
vere hysteresis upon temperature reversals near the phase 
boundary, is highly unusual. It is not correct to assume 
that simply not being in equilibrium would account for 
different critical behavior [37|. It is likely that such be- 
havior is more generic in systems undergoing phase tran- 
sitions in the presence of quenched disorder. Although 
FexZni_xF2 and its isomorphs in applied fields are the 
most characterized examples of the RFIM, other mag- 
netic 38, 39| and ferroelectric systems, as well as 
manganites |4l|, have been studied. 

SSRL is operated by Stanford University for the U.S. 
Department of Energy, Offlce of Basic Energy Sciences. 
The work at UCSC was funded by Department of Energy 
Grant No. DE-FG02-05ER46181. The work at Stanford 
was supported by the U.S. Department of Energy un- 
der Contract Nos. DE-FG03-99ER45773 and DE-AC03- 
76SF00515, by NSF Grant Nos. DMR-9400372 and 
DMR-9802737, and by the A.P. Sloan Foundation. We 
acknowledge Onuttom Narayan for useful discussions. 



12 



30 



20 - 



10 



a ^ A i 



H = ^lOT 
> lOT ^H=0 
□ H=0 ^lOT 
. lOT ^7.6T 
V 7.6T ^lOT 



_L 



_L 




4 6 
H(Tesla) 



10 



FIG. 18: Comparison of the field cycle Bragg scattering data 
at r = 63.10 K. Open symbols represent data from 'field- 
raising' and solid symbols represent data from 'field- lowering'. 
The dotted curves are guides to the eye. 
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